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SUMMARY 


Work  performed  under  the  subject  contract  (from  October  1977  to 
January  1980)  was  aimed  at  developing  the  methodology  for  computing 
unsteady  laminar  boundary  layers  with  flow  reversal  and  establishing  a 
new  unsteady  separation  criterion.  We  have  completed  calculation  of 
two-dimensional,  unsteady  boundary  layers  based  strictly  on  the  analogy 
with  the  steady,  three-dimensional  problems.  The  results  demonstrated 
the  validity  of  our  proposed  unsteady  separation  criterion.  As  a 
preclude  to  understanding  the  full,  three-dimensional,  unsteady  problem, 
we  also  treated  an  unsteady  symmetry-plane  boundary  layer.  Initial 
test  calculations  were  successful,  but  the  bulk  of  the  computation 
remains  to  be  carried  out. 
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I .  INTRODUCTION 


Unsteady  viscous  flow  phenomena  are  of  fundamental  importance  in 
many  practical  aerospace  applications.  Besides  occurring  from  laminar 
instability  and  transition  at  high  Reynolds  number,  unsteady  flow  may 
arise  either  when  the  free  stream  is  unsteady  or  when  the  body  experiences 
some  unsteady  motion.  Examples  of  the  latter  include  the  accelerating 
(or  decelerating)  motion  during  take-off  (or  landing)  operations  and 
the  periodic  motion  of  rotating  machineries  or  vehicles  (such  as 
helicopters,  turbines,  etc.). 

The  state  of  research  into  the  unsteady  aspects  of  the  laminar 
boundary  layer  theory  is  rather  unsatisfactory  compared  to  that  for  the 
steady  case.  Even  for  two-dimensional  problems,  complete  solutions  are 
still  scarce,  while  basic  questions  like  unsteady  separation  and  unsteady 
reversed  flow  remain  unsettled.  Erroneous  concepts  on  these  matters 
still  prevail  in  the  literature,  whereas  three-dimensional,  unsteady 
problems  are  virtually  untouched. 
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II.  RESEARCH  OBJECTIVES 


The  objective  of  this  research  program  was  to  investigate  the 
basic  aspects  of  two-dimensional,  and  then  three-dimensional,  unsteady 
viscous  flows.  To  obtain  concrete  results  and  definitive  conclusions 
without  excessive  efforts,  we  simplified  the  problems  by  following  the 
boundary  layer  approach  instead  of  the  semi-Navier-Stokes  or  full 
Navier-Stokes  approach.  We  particularly  focused  on  developing  a 
methodology  for  calculating  reversed  flow  and  establishing  a  workable 
criterion  for  determining  unsteady  separation. 

Although  the  problem  of  a  suitable  criterion  has  been  considered 
before  by  a  number  of  researchers,  we  felt  it  was  far  from  being  resolved. 
Consequently,  we  proposed  a  new  unsteady  separation  criterion  and 
devoted  most  of  our  efforts  during  the  contracting  period  to  demonstrating 
the  validity  of  the  proposed  criterion  with  specific  examples. 

Whereas  existing  separation  criteria  are  patterned  after  those 
for  two  dimensional,  steady  separation,  our  criterion  arises  from  a 
different  viewpoint:  Generalization  by  analogy  from  that  of  three- 
dimensional,  steady  separation. 

General  unsteady  separation  in  three-dimensions  was  also 
considered.  A  basic  requirement  for  the  general,  unsteady,  three- 
dimensinal  separation  criterion  is  that  it  must  consistently  Include 
both  the  known  steady,  three-dimensional  separation  criterion  and  the 
proposed  unsteady,  two-dimensional  separation  criterion  as  special 
cases.  Our  proposed  criterion  meets  these  requirements. 


III.  APPROACH 


We  used  the  analogy  approach  to  the  two-dimensional,  unsteady 
boundary  layer  problem,  i.e. ,  a  methodology  analogous  to  that  for 
solving  the  three-dimensional,  steady  boundary  layer  problem  was 
applied  to  solve  the  two-dimensional  unsteady  boundary  layer  problem. 

The  theoretical  basis  for  this  approach  was  the  observation  that  the 
respective  systems  of  governing  equations  embody  a  similar  mathematical 
structure.  Both  systems  contain  one  boundary-value  and  two  initial- 
value  independent  variables;  in  other  words,  they  are  hyperbolic  along 
two  coordinate  (including  time)  directions  and  elliptic  along  the  third 
(normal  to  the  body  surface)  direction.  Based  on  this  similarity,  we 
extended,  several  years  ago,  the  concept  of  the  zones  of  influence  and 
dependence  for  the  steady,  three-dimensional  case  to  the  unsteady,  two- 
dimensional  case,  even  though  the  number  of  equations  for  each  case  differed. 

In  the  present  work,  the  unsteady  analogy  was  further  explored 
on  several  fronts.  First,  the  comprehensive  computing  program  and 
associated  computational  logics  and  procedures  developed  previously  for 
the  steady,  three-dimensional  problems  were  converted  for  the  present, 
unsteady  two-dimensional  applications  by  simple  redesignations  of 
variables.  We  then  applied  the  concepts  for  calculating  reversed  cross- 
flow  over  an  inclined  body  of  revolution  (a  typical  steady,  three-dimensional 
problem)  to  calculation  of  the  unsteady,  two-dimensional  reversed  flow  — 
valid  as  long  as  the  respective  dependence  rules  are  satisfied.  An 
immediate  implication  of  the  calculabilitiy  of  unsteady  reversed  flow 
is  that  the  vanishing  of  skin  friction,  which  marks  the  beginning  of 
flow  reversal,  no  longer  signifies  separation  nor  automatically 
invalidates  the  boundary  layer  concept. 

The  unsteady  analogy  was  further  extended  to  deal  with  the  question  of 
unsteady  separation.  For  steady,  three-dimensional  flows,  separation 
line  has  been  defined  as  an  envelope  of  the  limiting  streamlines. 

This  envelope  criterion,  originally  based  on  flow  visualization  experiments, 
has  since  been  reaffirmed  by  available  numerical  solutions.  During 
the  present  research,  an  analogous  envelope  criterion  for  determining 
unsteady  separations  was  proposed  and  demonstrated. 
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Based  on  our  conviction  that  future  research  into  unsteady  flow 
should  focus  more  on  three-dimensional  problems,  we  also  undertook,  as 
a  first  step,  to  investigate  the  unsteady  boundary  layer  along  the 
symmetry  plane  of  an  inclined  body  of  revolution.  Previously,  we 
considered  the  steady  counterpart  and  obtained  results  that  revealed 
new  three-dimensional  features  of  fundamental  significance.  We  expect 
that  similar  situation  may  also  occur  for  the  unsteady  case.  Three- 
dimensional  unsteady  flow  characteristics  are  expected  to  be  significantly 
different  from  our  current  understanding  based  on  two-dimensional  or 
axisymmetrical  cases.  In  particular,  the  open  separation  phenomenon 
we  found  in  the  three-dimensional  steady  case  would  also  be  present  in 
the  corresponding  unsteady  case.  The  way  in  which  an  open  separation 
developed  as  the  boundary  layer  grows  from  an  impulsive  start  is,  for 
example,  very  intriguing  and  may  shed  some  light  on  the  previously 
unexplored  structure  of  three-dimensional  unsteady  flow. 
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IV.  DISCUSSION  AND  RESULTS 


From  our  investigations  of  two-dimensional  unsteady  separation, 
we  concluded  that  the  existing  unsteady  separation  criteria  (MRS 
criterion,  after  Moore,  Rott  and  Sears;  unsteady  Goldstein  singularity  (D 
criterion)  remain  unproven. 

Although  the  solutions  of  Williams  and  Johnson  (2)  demonstrated 
the  MRS  criterion,  the  special  assumptions  used  there  reduced  the 
problems  essentially  to  steady  ones  for  which  the  validity  of  the  MRS 
criterion  was  not  in  doubt.  The  cylinder  solution  by  Telionis  and 
Tsahalis  (3)  gave  support  to  the  Goldstein  singularity  criterion,  but 
this  solution  had  been  subsequently  contradicted  by  later  calculations. 

Apart  from  the  basic  validity  question,  both  of  these  criteria  are 
inconvenient  to  apply;  in  addition,  their  extension  to  three-dimensional 
cases  has  not  been  shown. 

To  analogously  apply  the  envelope  separation  criterion  for 
steady,  three-dimensional  flow  to  the  two-dimensional,  unsteady  flow, 
pertinent  equations  were  formally  transformed  from  one  to  the  other; 
similar  limiting  streamlines  were  defined  in  the  x,  t-plane  for  the 
unsteady  case;  and  the  envelope  of  those  limiting  streamlines  was 
proposed  as  a  new  unsteady  separation  criterion.  This  criterion 
determines  separation  without  relying  on  usual  indications,  such  as 
rapid  increase  of  boundary  layer  thickness,  etc.,  and  it  is  also 
consistent  with  the  criterion  for  general  unsteady  separation  in  three 
dimensions. 

Analogous  to  the  three-dimensional,  steady  case,  separation  for  the 
two-dimensional,  unsteady  case  can  also  be  classified  as  "open"  or  "closed”, 
i.e.,  the  separation  line  defined  as  an  envelope  of  the  limiting 
streamlines  does  not  or  does  divide  the  x,  t-plane  into  two  unconnected 
regions.  To  illustrate  these  types  of  separation,  two  previously 
studied  examples  were  recalculated:  unsteady  response  to  a  sudden 
change  of  the  Howarth  steady  problem  (closed  separation)  and  unsteady 
growth  for  an  impulsively  started  circular  cylinder  from  rest  (open  separation). 
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Results  confirmed  envelope  formation  by  analogous  unsteady 
limiting  stream  lines,  and  the  separation  so  determined  was  compared  to 
that  of  Telionis  and  Tsahalis'.  For  the  first  problem,  (closed 
separation)  agreements  were  noted  in  general  trends  of  the  separation 
line  in  relation  to  the  zero  skin-friction  line,  but  discrepancies  were 
found  in  actual  locations  of  these  lines.  Also  the  Telionls-Tsahalis 
results  showed  a  much  faster  approach  toward  the  final  steady-state 
condition  than  did  ours.  For  the  second  problem  (open  separation), 
results  for  the  zero  skin-friction  line  agreed  very  well,  but  those  for 
the  separation  line  were  entirely  different.  Our  separation  was  determined 
by  convergence  to  an  envelope  of  the  limiting  streamlines,  theirs  by 
the  singularity  criterion.  Our  calculations  detected  no  singularity  at 
the  place  and  time  shown  by  their  calculations.  Moreover,  we  predicted 
separation  at  larger  times. 

Our  prediction  appeared  to  be  in  good  agreement  with  that  of  Van 
Dommelen  and  Shen  based  on  a  Lagrangian  approach.  Although  our  results 
for  the  vanishing  skin  friction  point  and  the  boundary  layer  thickness 
also  agreed  well  with  those  by  Cebeci  (5)^  yet  the  separation  conclusions 
reached  were  very  different.  He  contended  that  no  separation  occurred 
at  any  finite  time,  we  concluded  just  the  opposite.  Interestingly, 
however,  both  we  and  Cebeci  claimed  to  agree  with  the  general  conclusion 
of  Proudman  and  Johnson  (6). 

We  also  described  general  unsteady  separation  in  three  dimensions 
based  on  the  application  of  a  steady,  three-dimensional  criterion  at 
successive  instants  of  time.  This  method  simulates  the  usual  sequence 
of  observing  an  unsteady  flow,  and  does  not  necessarily  imply  a  quasi¬ 
steady  approach.  Alternatively,  it  may  be  seen  as  applying  the  unsteady, 
two-dimensional  criterion  to  staggering  fixed  coordinate  planes.  Hence, 
the  three-dimensional,  unsteady  separation  we  proposed  includes  both 
the  steady,  three-dimensional  separation  and  the  unsteady,  two-dimensional 
separation  as  special  cases.  Our  unsteady,  three-dimensional  separation 
criterion  differed  from  that  of  Eichelbrenner  (7)  which  contained  the  closed 
type  of  separation  but  not  the  open  type. 


The  details  of  our  work  were  documented  in  the  report: 

Wang,  K.C.,  "Unsteady  Boundary  Layer  Separation,”  Martin  Marietta 
Laboratories  TR  79-16c,  April  1979. 

Subsequent  calculated  results  and  comparisons  with  other  work  will 
soon  be  incorporated  into  a  revised  version  for  journal  publication. 

Work  on  the  unsteady  boundary  layer  along  the  symmetry  plane  of 
an  inclined  ellipsoid  of  revolution  is  well  under  way.  The  detailed 
formulation,  change  of  computing  program,  and  generation  of  initial 
profiles  have  all  been  completed.  Initial  calculations  have  been 
successful,  and  different  grid  steps  in  time  as  well  as  along  the 
symmetry  plane  have  been  tested  to  determine  optimal  mesh  sizes  for 
conserving  computing  time.  However,  the  bulk  of  the  calculations  for 
this  work  remains  to  be  carried  out.  Table  I  illustrates  the  recently 
calculated  skin  frictions  on  both  the  windside  and  the  leeside  of  the 
body  at  small  times. 
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TABLE  I.  SKIN  FRICTIONS 
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